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I.  INTRODUCTION 


One  of  the  most  important  aerodynamic  performance  characteristics  for 
shell  is  the  total  drag.  The  total  drag  for  projectiles  can  be  divided  into 
three  components:  (1)  pressure  drag  (excluding  the  base),  (2)  viscous  (skin 

friction)  drag,  and  (3)  base  drag.  The  base  drag  is  a  major  contributor  to 

the  total  drag  and  can  be  as  much  as  50-60%  of  the  total  drag.  It  is,  thus, 

important  to  minimize  this  part  of  the  drag. 

An  effective  means  of  reducing  the  base  drag  is  by  base  bleed.  In  this 
method,  a  relatively  small  quantity  of  low  velocity  gas  is  injected  into  the 
dead  air  region  immediately  behind  the  base.  This  changes  the  structure  of 
the  flow  field  in  the  base  region  (See  Figure  1)  which  results  in  an  increase 
in  the  base  pressure  or  reduction  in  the  base  drag.  As  mass  flow  rate  of  gas 
injected  into  the  base  region  is  increased  from  zero,  the  base  pressure  first 
increases  (See  Figure  2)  until  a  maximum  is  reached.  With  further  increase  in 
mass  injection,  the  base  pressure  falls  to  a  minimum  and  then  rises  again. 
The  'base  bleed'  region  is  the  first  part  of  this  curve  up  to  the  maximum 
which  corresponds  to  very  small  rates  of  gas  injection  (of  the  order  of  a  few 
percent). 

The  drag  reduction  due  to  base  bleed  at  supersonic  speeds  is  of  practical 
importance.  The  concept  of  mass  injection  at  the  base  has  received  consid¬ 
erable  interest  in  the  past  and  an  excellent  review  is  reported  by  Murthy  et 
al.1  As  reported  in  Reference  1,  there  is  a  lack  of  detailed  experimental 
measurements  of  the  base  flow  phenomena  with  or  without  mass  injection.  Super¬ 
sonic  flow  over  cylindrical  afterbodies  with  gas  injection  has  been  experi¬ 
mentally  studied  by  Bowman  and  Clayden.2  Their  results  showed  the  effect  of 
mass  injection  on  the  base  pressure  for  a  range  of  Mach  numbers  between  1.5 
and  3.0.  For  low  Mach  numbers  (M  «  1.5),  an  increase  in  base  pressure  was 
obtained  with  increasing  mass  flow  rate  of  gas  injection.  For  higher  Mach 
numbers,  the  typical  pattern  of  the  effect  of  gas  injection  on  the  base  pres¬ 
sure  shown  in  Figure  2  was  obtained.  The  effect  of  heating  the  injected  gas 
results  in  a  further  but  relatively  small  increase  in  base  pressure  or  de¬ 
crease  in  base  drag.3  More  recently.  Schilling4  performed  an  experimental 
study  to  investigate  base  bleed  effect  on  cylindrical  and  boattailed  after¬ 
bodies  where  primary  emphasis  was  to  assess  the  effect  of  the  tail  fins  on  the 


1.  Murthy,  S.N.B. ,  (Ed.),  Progress  in  Astronautics  and  Aeronautics:  Aero¬ 
dynamics  of  Base  Combustion,"  Vol .  40,  AIAA,  Dew  York,  1976. 

2.  Bowman,  J ,E .  and  Clayden,  W.A.,  "Cylindrical  Afterbodies  in  Supersonic 
Flow  with  Gas  Ejection,"  AIAA  Journal,  Vol.  5,  No.  8,  August  1967,  pp. 
1524-1525. 

5.  Clayden,  W.A.  and  Bowman,  J.E.,  "Cylindrical  Afterbodies  at  =  2  with 
Hot  Gas  Ejection,"  AIAA  Journal,  Vol.  6,  No.  12,  December  1968,  pp.  2429- 
2431  . 

4.  Schilling,  H.  "Experimental  Investigation  on  the  Base-Bleed-Effect  for 

Body-Tail-Combinations,"  Proceedings  of  the  8th  International  Symposium  on 
Ballistics,  Amsterdam,  Holland,  1984. 


effectiveness  of  base  bleed.  These  experiments  showed  a  strong  interaction 
between  the  fins  and  the  base  region  flow  and  the  negative  effect  it  has  on 
the  base  drag  reduction.  Questions  regarding  the  difference  between  hot  gas 
injection  and  cold  gas  injection,  the  effect  of  base  bleed  on  stability  char¬ 
acteristics,  the  effect  of  spin  on  the  base  bleed  effect  and  base  combustion 
still  remain  to  be  answered.  Very  few  semi -empi ri cal  methods  are  available  to 
predict  the  reduction  in  base  drag  due  to  base  bleed.  One  such  method  has 
been  developed  by  Hellgren5  in  Sweden.  An  attempt  has  been  made  in  this 
method  to  take  into  account  the  variation  of  burning  conditions  in  the  base 
bleed  motor.  Uncertainty  in  the  base  drag  coefficient  and  also  in  the  effect 
of  spin  on  the  burning  base  bleed  propellant  exists  in  this  method.  The 
application  of  sophisticated  computational  techniques  to  the  practical  base 
bleed  problem  is  in  its  infancy.  Limited  computational  work  has  been  reported 
by  Sul  1 i ns  et  al6  which  dealt  with  the  numerical  computation  of  base  region 
flow  with  parallel  gas  injection  using  two-dimensional  Navier-Stokes  equa¬ 
tions.  Recently,  a  new  numerical  capability  has  been  developed  by  Sahu  et  al7 
to  compute  the  base  region  flow  using  Azimuthal-Invariant  thin-layer  Navier- 
Stokes  equations.  Since  the  entire  projectile  flow  field  was  computed,  it  was 
possible  to  obtain  the  total  drag  at  transonic  speeds.  An  attempt  was  also 
made  to  include  the  base  injection  at  transonic  speeds  using  crude  base  bound¬ 
ary  conditions  although  the  primary  interest  of  the  base  bleed  concept  is  in 
the  supersonic  speed  regime. 

This  paper  describes  a  computational  investigation  of  the  effect  of  cen¬ 
tered  base  bleed  on  the  base  region  flow  field  and  on  the  base  pressure  at 
supersonic  speeds  using  basically  the  same  numerical  procedure  of  Reference  7. 
A  new  interactive  base  bleed  boundary  condition  procedure  has  been  developed 
to  determine  the  bleed  exit  boundary  conditions.  The  resulting  numerical 
capability  is  used  to  compute  supersonic  flow  over  two  axi symmetric  cylindric¬ 
al  afterbodies  with  gas  injection  and  the  computed  results  are  compared  to 
experiment.  This  is  a  first  step  towards  obtaining  a  capability  to  predict 
the  effect  of  base  bleed  (including  combustion  in  the  base  region)  on  the 
aerodynamics  of  Army  projectiles. 


5.  Hellgren,  R.V.,  "Range  Calculation  for  Base  Bleed  Propellants,  "  Proceed¬ 
ings  of  the  6th  International  Symposium  on  Ballistics,  Orlando,  Florida, 
1981 . 

6.  Sullins,  G.A.,  Anderson,  J.D.,  and  Drummond,  J.P.,  "Numerical  Investiga¬ 
tion  of  Supersonic  Base  Flow  with  Parallel  Injection,"  ,41,4,4  Paper  No. 
82-1001,  June  1982. 

7.  Sahu,  J.,  Nietubicz,  C.J.,  and  Steger,  J.L.,  ,rNavier-Stokes  Computations 
of  Projectile  Base  Flow  with  and  without  Base  Injection,"  US  Army 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  ARBRL- 
TR-02532,  November  1983.  (AD  A135738)  (Also  see  AIAA  Journal,  Vol.  23, 
No.  9,  September  1985,  pp.  1348-1355 ) 


II.  COMPUTATIONAL  TECHNIQUE 


The  Azimuthal  Invariant  (or  Generalized  Axi symmetri c)  thin-layer  Navier 
Stokes  equations  for  general  spatial  coordinates  £,  n,  c  can  be  written 
as:8 


aTq  +  a^E  +  a;G  +  h  =  Re_1a;$ 


(1 


where  £  =  £(x,y,z,t)  is  the  longitudinal  coordinate 

n  =  n(y,z,t)  is  the  circumferential  coordinate 
C  =  C(x,y,z,t)  is  the  near  normal  coordinate 
t  =  t  is  the  time 
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8.  Nietubicz ,  C.J.,  Pulliam,  T.H.,  and  Steger,  J.L.,  "Numerical  Solution  of 
the  Azimuthal-Invariant  Navier-Stokes  Equations, "  US  Army  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  ARBLR-TR-02227 , 
March  1980.  (AD  A0857161  (Also  see  AIAA  Journal,  Vol.  18,  No.  12, 
December  1980,  pp.  1411-1412 )  - 
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"<{x  *  'y  +  Ez>uc  *  (p/3){cxuc  *  Cyv{  * 
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+  (p/3) (cxu  +  cyv  +  ^zw)(4xy^  +  cyv?  +  czwc)} 


The  velocities 


u  =  h  +  5xu  +  V  +  CzW 

v  =  nt  +  oxu  +  nyv  +  n2w  (2) 

w  -  «t  +  y  +  y  +  y 

represent  the  contravari ant  velocity  components. 

The  Cartesian  velocity  components  (u,  v,  w)  are  nondimensional i zed  with 
respect  to  am  (free  stream  speed  of  sound).  The  density  (p)  is  referenced 
to  pa  and  total  energy  (e)  to  p^a*.  The  local  pressure  is  determined  using 
the  equation  of  state, 

P  =  (y  -  l)[e  -  0.5p(u2  +  V2  +  w2)]  (3) 

where  y  is  the  ratio  of  specific  heats. 

In  Equation  (1),  axisymmetric  flow  assumptions  have  been  made  which  re- 

A 

suit  in  the  source  term,  H.  The  details  of  how  this  is  obtained  can  be  found 
in  Reference  (8)  and  are  not  discussed  here.  Equation  (1)  contains  only  two 
spatial  derivatives.  However,  it  retains  all  three  momentum  equations  and 
allows  a  degree  of  generality  over  the  standard  axisymmetric  equations.  In 
particular,  the  circumferential  velocity  is  not  assumed  to  be  zero  thus  allow¬ 
ing  computations  for  spinning  projectiles  to  be  accomplished.  This  is 
especially  important  to  study  the  effect  of  spin  on  the  base  bleed  effect. 

The  numerical  algorithm  used  is  the  Beam-Warming  fully  implicit,  approxi¬ 
mately  factored  finite  difference  scheme.  The  algorithm  can  be  first  or 


wr 


second  order  accurate  in  time  and  second  or  fourth  order  accurate  in  space. 
Since  the  interest  is  only  in  the  steady-state  solution.  Equation  (1)  is 
solved  in  an  asymptotic  fashion  and  first  order  accurate  time  differencing  is 
used.  The  spatial  accuracy  is  fourth  order.  Details  of  the  algorithm  are 
included  in  References  9-11. 

To  suppress  high  frequency  components  that  appear  in  regions  containing 
severe  pressure  gradients  e.g.,  shocks  or  stagnation  points,  artificial  dis¬ 
sipation  terms  are  added.  In  the  present  application,  a  switching  dissipation 
model  is  used  which  is  a  blend  of  second  and  fourth  order  dissipation  terms. 
This  is  similiar  to  the  model  used  by  Pulliam12  and  uses  a  fourth  order  dis¬ 
sipation  in  smooth  regions  and  switches  to  a  second  order  dissipation  in 
regions  containing  high  pressure  or  density  gradients.  Incorporation  of  this 
dissipation  model  has  resulted  in  the  improvement  of  the  quality  of  the 
results  and  has  made  the  code  more  robust. 

The  numerical  code  used  computes  the  full  flow  field  over  a  projectile  or 
a  missile  including  the  base  region  using  a  unique  flow  field  segmentation 
procedure.  An  important  advantage  of  this  procedure  lies  in  the  preservation 
of  the  sharp  corner  at  the  base.  The  details  of  these  can  be  found  in  Refer¬ 
ence  7.  For  the  computation  of  turbulent  base  flows,  the  two-layer  algebraic 
Bal dwi n-Lomax  turbulence  model 13  is  used.  Higher  order  two-equation  or  more 
sophisticated  turbulence  models  need  to  be  considered  for  such  flows  in 
future. 


III.  BOUNDARY  CONDITIONS 

The  no  slip  boundary  condition  for  viscous  flow  is  enforced  by  setting 
the  contravariant  velocities  to  zero  i.e. ,  U  =  V  =  W  =  0  on  the  body  surface. 
Inviscid  boundary  conditions  are  used  at  the  base  for  the  case  of  no  gas 
injection.  The  pressure  on  the  body  surface  is  obtained  by  solving  a  combined 
momentum  equation  and  density  is  extrapolated.  First  order  extrapolation  is 
used  for  both  pressure  and  density  at  the  base.  Along  the  computational  cut, 
the  flow  variables  above  and  below  are  averaged  to  determine  the  conditions  on 


9.  Steger,  J.L.,  "Implicit  Finite  Difference  Simulation  of  Flow  About 
Arbitrary  Geometries  with  Application  to  Airfoils,"  AIAA  Journal,  Vol. 

16,  No.  4,  July  1978,  pp.  679-686. 

10.  Pulliam,  T.H.  and  Steger,  J.L.,  "On  Implicit  Finite-Difference  Simula¬ 
tions  of  Three-Dimensional  Flow,  "  AIAA  Journal,  Vol .  18,  No.  2,  February 
1980,  pp.  159-167. 

11.  Beam,  R.  and  Warming,  R.F.,  "An  Implicit  Factored  Scheme  for  the  Compres¬ 
sible  Navier-Stokes  Equations,"  AIAA  Paper  No.  77-645,  June  1977. 

12.  Pulliam,  T.H.,  "Artificial  Dissipation  Models  for  the  Euler  Equations," 
AIAA  Paper  No.  85-0438,  January  1985. 

13.  Baldwin,  B.S.  and  Lomax,  H.,  "Thin-Layer  Approximation  and  Algebraic 
Model  for  Separated  Turbulent  Flows,"  AIAA  Paper  No.  78-257,  1978. 


the  cut.  A  symmetry  boundary  condition  is  imposed  on  the  centerline  of  the 
wake  region  with  first  order  extrapolation  for  density  and  pressure.  Free 
stream  conditions  are  imposed  on  the  outer  boundary  while  extrapolated  values 
of  the  variables  are  used  at  the  downstream  boundary. 

For  the  base  bleed  case,  boundary  conditions  along  the  base  where  gas  is 
injected  need  to  be  modified.  Here,  a  small  amount  of  mass  is  injected  into 
the  near  wake  behind  the  center  part  of  the  base.  The  amount  of  gas  injection 
is  usually  defined  in  terms  of  a  mass  injection  parameter. 


I 


m 


where  m.  is  the  mass  flow  at  the  bleed  exit  and  Ab  is  the  area  at  the  base. 

•J 

At  the  bleed  exit,  the  flow  is  subsonic  and  three  quantities  need  to  be  speci¬ 
fied.  Total  temperature  is  used  at  the  exit  along  with  an  extrapolated  value 
of  static  pressure.  The  static  pressure  at  the  exit  is  set  equal  to  local 
pressure  which  is  a  characteristic  of  subsonic  flows.  Such  extrapolated 
static  pressure  boundary  conditions  have  been  used  in  other  numerical  computa¬ 
tions.1^  15  In  addition,  the  value  of  total  pressure  at  the  exit  is  needed  to 
define  all  the  boundary  conditions.  This  is  determined  iteratively  to  satisfy 
the  required  mass  injection  parameter  in  the  following  manner.  Knowing  the 
total  temperature  (Toj)  and  mass  injection  parameter  (I), 

(a)  assume  total  pressure,  POj 

(b)  obtain  Pj  by  extrapolation 

(c)  calculate  Mach  number,  M,-  using  Po./P.-  =  (1  +  .2M.2)3'5 

J  J  J  J 

(d)  compute  fj  using  Toj/Tj  =  1  +  . 2Mj 2 

(e)  obtain  the  velocity  u .=  M.  a.  =  M-  /yRT . 

J  J  J  J  J 

(f)  obtain  the  density  Pj  from  the  equation  state,  Pj  =  PjRTj 

(g)  compute  I  =  p  j  u  j  A  j /P®u«»Ab 


14.  Shrewsbury,  G.D.,  "Analysis  of  Circulation  Control  Airfoils  Using  an 
Implicit  Navier-Stokes  Solver,"  AIAA  Paper  NO.  85-0171,  January  1985. 


15. 


Drummond,  J.P.,  ’Numerical  Study  of  a  Ramjet  Dump  Combustor  Flow  Field,  " 
AIAA  Paper  No.  83-0421,  January  1983. 


If  the  value  of  this  I  is  not  the  desired  value,  then  we  go  back  to  step  (a) 
and  continue  the  iteration  process.  This  iteration  procedure  is  repeated  at 
each  time  step.  The  total  pressure  and  the  total  temperature  are  assumed  to 
be  constant  across  the  bleed  exit. 


IV.  RESULTS 

Computations  of  the  flow  field  for  two  cylindrical  afterbodies  with  base 
bleed  have  been  made  at  supersonic  speed  and  at  zero  angle  of  attack.  The 
first  set  of  computed  results  correspond  in  part  to  the  experimental  investi¬ 
gation  on  the  base  bleed  effect  by  Schilling4  whereas  the  second  set  of 
results  are  compared  to  the  experimental  measurements  by  Bowman  et  al 2  at  the 
Royal  Armament  and  Development  Establishment  (RARDE),  England. 

The  experimental  investigation  by  Schi 1 1 ing4  was  conducted  at  the  Tri- 
sonic  Wind  Tunnel  at  the  Deutsche  Forschungs-  und  Versuchsanstalt  fur  Luft- 
und  Raumfahrt  e.  V.  Kdln  (DFVLR)  in  West  Germany  for  a  missile  configuration. 
A  computational  grid  for  this  case  showing  the  cylindrical  afterbody  region  is 
shown  in  Figure  3.  The  full  grid  consisted  of  156  points  in  the  streamwise 
direction  and  60  points  in  the  normal  direction.  The  grid  has  been  adapted  to 
the  free  shear  layer  as  the  solutions  developed.  The  same  grid  was  used  for 
calculations  with  and  without  base  bleed.  The  minimum  grid  spacing  near  the 
wall  and  in  the  free  shear  layer  is  .00002  caliber.  Numerical  computations 
were  made  at  M^,  =  1.7  and  2.5  at  zero  angle  of  attack. 

Figure  4  shows  the  velocity  vectors  in  the  base  region  for  the  case  of  no 
gas  injection  whereas  Figures  5  and  6  show  the  velocity  vectors  in  the  base 
region  with  gas  injection  for  Mw  =  2.5.  In  Figure  4,  the  reci rcul atory  flow 

field  in  the  near  wake  is  clearly  evident.  The  flow  expands  at  the  base  and 
the  reattachment  point  is  about  .75  caliber  downstream  of  the  base.  The  ef¬ 
fect  of  gas  injection  at  the  base  is  shown  in  Figures  5  and  6  for  mass  injec¬ 
tion  rates  I  =  .01  and  .02,  respectively.  As  shown  in  Figure  5  for  I  =  .01, 
the  separation  bubble  (for  the  no  bleed  case)  has  been  displaced  downstream 
and  the  size  of  the  separation  bubble  has  also  been  reduced.  The  flow  field 
in  the  near  wake  has  changed  considerably.  There  are  two  stagnation  points, 
one  at  X/D  -  15  and  the  other  at  X/D  =  15.4.  This  is  a  typical  flow  field 
pattern  with  base  bleed  at  a  low  injection  rate.  As  I  is  increased  to  .02, 
more  mass  is  injected  into  the  near  wake  and  strongly  affects  the  base  region 
flow.  As  seen  in  Figure  6,  the  separation  bubble  that  was  displaced  down¬ 
stream  somewhat  at  I  =  .01  has  now  been  eliminated  and  the  beginning  of  a 
small  separation  region  near  the  base  is  seen.  For  this  injection  rate,  the 
Mach  number  at  the  bleed  exit  approaches  sonic  value.  Increasing  the  mass 
injection  rate  further  gives  Mach  number  unity  at  the  bleed  exit  and  results 
in  the  separation  region  in  the  near  wake  which  resembles  that  of  a  jet 
effect. 

Figure  7  shows  the  pressure  distribution  as  a  function  of  the  longitudi¬ 
nal  position  over  the  cylindrical  afterbody  and  along  the  upper  shear  layer  in 
the  wake  region.  The  strong  expansion  at  the  base  corner  and  the  recompres¬ 
sion  waves  downstream  of  the  base  corner  are  clearly  seen  for  the  zero  base 
bleed  case  (I  =  0).  The  pressure  distribution  for  the  base  bleed  case  is 
presented  in  this  figure  for  I  =  .01  and  I  =  .02.  As  I  increased  to  .01  from 
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zero,  the  expansions  at  the  base  and  the  recompression  downstream  of  the  base 
are  both  weakened.  Increasing  I  further  has  the  effect  of  weakening  the  waves 
even  more. 

The  effect  of  mass  injection  on  the  base  pressure  and  hence  on  base  drag 
is  of  primary  interest.  Figure  8  shows  the  base  drag  as  a  function  of  the 
mass  injection  parameter  for  M^  =  1.7  and  a  =  0.  1=0  corresponds  to  the  case 

of  no  base  bleed.  As  I  is  increased,  the  base  pressure  increases  in  this  case 
and  thus,  base  drag  is  reduced.  This  is  true  for  both  I  =  .01  and  I  =  .02. 
The  computed  base  drag  is  compared  to  the  experimental  data  and  as  can  be 
seen,  the  computed  results  are  overpredicted  by  about  15%.  The  trend  of  base 
drag  reduction  seen  experimentally  however,  has  been  clearly  predicted  by  the 
computation.  Additionally,  the  Mach  number  at  the  bleed  exit  is  less  than 
unity  for  all  the  mass  injection  rates  considered  here.  The  results  for  = 

2.5  is  shown  in  Figure  9.  Here  the  computed  base  drag  is  compared  to  experi¬ 
ment  and  excellent  agreement  is  obtained  for  the  zero  bleed  condition  (I  =  0). 
As  the  mass  flow  rate  is  increased  to  .01,  base  pressure  increases  and  base 
drag  is  reduced.  Comparison  between  the  computation  and  the  experiment  again 
indicates  10-15%  disagreement.  The  Mach  number  at  the  bleed  exit  (Mj)  is 

approaching  unity.  As  I  is  increased  further  to  .02,  Mj  has  reached  the  sonic 

condition  and  stays  sonic  as  I  is  increased  further.  Here  the  trend  reverses 
and  a  small  increase  in  base  drag  is  found.  For  I  =  .025,  the  computed  base 
drag  is  in  good  agreement  with  the  experiment. 

Supersonic  flow  for  another  cylindrical  afterbody  with  base  bleed  has 
been  computed.  Experimental  measurements  for  this  case  have  been  made  by  Bow¬ 
man  et  al2  at  RARDE,  England.  The  model  configuration  was  a  projectile  con¬ 
sisting  of  a  two  caliber  ogive  nose  and  a  five  caliber  cylindrical  afterbody. 
The  bleed  exit  diameter  was  equal  to  0.4  base  diameter.  Computations  for  this 
case  have  been  made  at  two  supersonic  speeds.  M^  =  1.88  and  M^  =  2.48  at  zero 

angle  of  attack.  The  computational  grid  used  in  this  case  for  =  1.88  is 

shown  in  Figure  10.  The  grid  is  adapted  to  the  free  shear  layer  in  the  wake 
and  clustering  of  grid  points  is  made  near  the  base  bleed  exit.  The  number  of 
grid  points  in  the  streamwise  direction  is  130  and  40  points  are  used  in  the 
normal  direction.  The  qualitative  results  obtained  in  the  base  region  with 
base  bleed  are  similar  to  that  discussed  earlier.  The  effect  of  gas  injection 
on  the  base  pressure  is  discussed  next. 

Figures  11  and  12  show  the  ratio  of  base  pressure  to  free  stream  static 
pressure  as  a  function  of  the  mass  injection  parameter,  I  for  M^  =  1.88  and 

2.48,  respectively.  As  shown  in  Figure  11  for  M^  =  1.88,  base  pressure  in¬ 
creases  as  I  is  increased  from  zero.  The  agreement  between  the  computed 
result  and  the  experiment  is  very  good.  As  I  =  .02  is  approached,  the  base 
pressure  ceases  to  rise  according  to  the  experimental  data.  Figure  12  shows 
the  base  pressure  with  and  without  base  bleed  between  the  computation  and  the 
experiment.  As  I  is  increased  from  zero  to  .01,  both  experiment  and  the  comp¬ 
utation  show  an  increase  in  the  base  pressure  and  reasonably  good  agreement  is 
obtained.  Further  increase  in  I  clearly  shows  a  large  decrease  in  the  base 
pressure  as  is  seen  in  the  experimental  data.  The  computation  doesn't  show 
that  much  of  a  decrease  in  the  base  pressure  although  the  trend  is  in 


agreement  with  the  experimental  data.  Also,  the  flow  was  assumed  to  be  sub¬ 
sonic  at  the  bleed  exit  in  the  computation;  however,  it  is  not  clear  what  the 
exit  condition  was  in  the  experiment.  The  comparison  of  the  results  for  this 
case  must  be  viewed  in  that  light. 

V.  CONCLUDING  REMARKS 

A  computational  study  has  been  made  for  base  region  flow  over  two  cylin¬ 
drical  afterbodies  with  base  bleed  at  supersonic  speeds.  An  iterative  boun¬ 
dary  condition  procedure  was  developed  for  the  base  bleed  effect  and  was  used 
to  show  the  effect  of  base  bleed  on  base  pressure  or  base  drag.  The  thin- 
layer  form  of  the  compressible  Navi er-Stokes  equations  was  solved  using  a 
time-dependent  implicit  numerical  algorithm.  Numerical  results  show  the 

qualitative  effect  of  the  mass  injection  on  the  near  wake  flow  field.  The 
expansions  at  the  base  and  the  recompression  downstream  of  the  base  are 

weakened  by  the  gas  injection.  For  low  supersonic  Mach  number  -  1.75,  base 
pressure  rises  smoothly  as  injected  mass  flow  rate  is  increased.  For  the 

higher  Mach  number  =  2.5,  increase  in  base  pressure  or  decrease  in  base  drag 
is  obtained  at  lower  mass  flow  rates  (I  =  .01).  Further  computational  invest¬ 
igation  is  needed  to  include  the  effect  of  hot  gas  injection  at  the  base. 
Future  work  will  be  directed  towards  obtaining  a  predictive  capability  to  com¬ 
pute  projectile  base  region  flow  with  combustion  for  practical  applications. 
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Figure  1.  Schematic  Illustration  of  Base  Flow  with  Base  Bleed 
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Figure  2.  Base  Pressure  vs  Mass  Flow  Rate  of  Injection 
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Figure  10.  Adapted  Grid  in  the  Base  Region  for  RARDE  Case,  =  1.88,  a 
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